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У статті, за допомогою розробленого нового інтегрального перетворення, знайдено температурне 
поле порожнього кусково-однорідного кругового циліндра в напрямку полярного радіуса, який 
обертається з постійною кутовою швидкістю навколо осі 07, з урахуванням кінцевої швидкості 
поширення тепла, у вигляді збіжних ортогональних рядів по функціях Бесселя і Фур'є. 
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І ре агіїсіе Ше деуеїоред пем/ іпіеєта! (гап5їогттайоп Їопйпа Ше (еппрегаїиге Пед егаріу ріесемізе 
ротобепеоц5 сігсиціаг суппаєг (омагаз8 Фе роїаг гадїця Па! гогагез аї а соп5іапі апешіаг угіосіїу ахіз8 О7, (аКіпе, 
іпіо ассоцпі їре бпіїе 5реед ої ргорабайоп ої реа іп (ре Їогт ої сопуегєепі огібоєгопаї зегіез іп Вез5е! 
Пшсбопя апа Боигіег 15. 

Кеу ууогаз: Пе іпіезта! Гаріасе ігап5їогтя, Еойигіег, геіахайоп те. 


Постановка проблеми, аналіз останніх досліджень і публікацій 

У феноменологічній теорії теплопровідності передбачається, що швидкість 
поширення тепла є нескінченно великою |1,2|. Однак при високих інтенсивних 
нестаціонарних процесах, що спостерігаються, наприклад, при вибухах, надзвукових 
потоках, великих швидкостях обертання, вплив скінченності величини швидкості 
поширення тепла на теплообмін стає помітним | 1,2|. 

В П| отримано узагальнене рівняння переносу енергії для рушійного елемента 
суцільного середовища, з урахуванням скінченності величини швидкості поширення тепла. 

Метою роботи є розробка нової узагальненої математичної моделі розподілу 
температурних полів у порожньому кусково-однорідному циліндрі у вигляді крайової 
задачі математичної фізики для рівняння теплопровідності, та розв'язання отриманої 
крайової задачі, розв'язок якої використовується підчас керування температурними полями. 

Як показує огляд літератури, теплообмін у циліндрах, які обертаються, вивчений 
на даний час ще недостатньо |1|. Показано, що чисельні методи дослідження 
нестаціонарних неосесиметричних задач теплообміну циліндрів, які обертаються, є не 
завжди ефективними, якщо мова йде про обчислення при великих швидкостях 
обертання. Так, доводиться |1|, що умови стійкості обчислень у методі кінцевих 
елементів і методі кінцевих різниць, що застосовуються до розрахунку нестаціонарних 
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неосесиметричних температурних полів циліндрів, які обертаються, визначаються 
аналогічними характеристиками. Ці умови мають вигляд: 


Розд «бо УВА. 


де Бу - критерій Фур'є, Ра - критерій Предводітелева. 

Якщо Ра« 10? , що відповідає кутовій швидкості обертання металевого циліндра 
оа1б71сек"! радіусом 100 мм, то змінні Дф та ДЕ, повинні бути підпорядковані 
таким умовам: 

Дра24107 ї ДЕ, «2-10719, 

Для рівномірно охолоджуваного циліндра, за умови Ві - 5 (Ві - критерій Біо), час, 
необхідний для того, щоб температура досягла 9095 стаціонарного стану, дорівнює 
Ко. 0.025(1). Це означає, що потрібно принаймні здійснити 1.3-10З операцій по часу 
для того, щоб було досягнуто стаціонарний розподіл температури. 

Більш того, потрібно відзначити, що протягом одного циклу обчислень потрібно 
здійснити 3.14-10? обчислень, так як внутрішній стан у кільці характеризується 


3.14-10? точками. У результаті видно, що число обчислень, необхідних для отримання 
чисельного результату, видається нереальним. 

Тому для вирішення крайової задачі, яка виникає при математичному 
моделюванні нестаціонарних процесів теплообміну в циліндрах, які обертаються, 
будемо застосовувати інтегральні перетворення. 

Постановка задачі 

Розглянемо розрахунок нестаціонарного температурного поля порожнього 
кругового циліндра зовнішнього радіуса К в циліндричній системі координат (», Ф, .), 
кусково-однорідного в напрямку полярного радіуса г, який обертається з постійною 
кутовою швидкістю (? навколо осі О7, з урахуванням кінцевої швидкості поширення 
тепла. І теплофізичні властивості якого в кожному шарі не залежать від температури за 
умови ідеального теплового контакту між шарами, а внутрішні джерела тепла відсутні. 
У початковий момент часу температура циліндра постійна Сі, а на зовнішній і 


внутрішній поверхні циліндра температура відома і не залежить від часу Сі С (ф) 
відповідно. 

Відносну температуру циліндра д(р, Ф, ') можна представити у вигляді: 
д(р,рі) якщо ре (ро,рі) 


(1) 
8,(р,р.і) якщо ре (рі, р») 


Фонрьг)- 
Відносні температури д.(р.д.1) 5-го шара циліндра | обчислюються за 
формулами: 


Т.Г о, Ф31)- С 
9 (р.фи)- Мо 0. 0 
тах 


9 
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де Т(р,т.г) - температури 5-го шару циліндра; Тиах- Максимальна температура 


циліндра; р ке ; 51,2. 


Розв'язок задачі 

В П| отримано узагальнене рівняння переносу енергії для рушійного елемента 
суцільного середовища, з урахуванням скінченності величини швидкості поширення 
тепла. Згідно з |1|, узагальнене рівняння балансу енергії твердого тіла, яке обертається 
з постійною кутовою швидкістю и»? навколо осі О/, теплофізичні властивості якого не 
залежать від температури, а внутрішні джерела тепла відсутні набуває вигляду: 


ДИ ЗД Ко З с ср.1ОГ. З0ч От 
усі ФО НТ, оо | зо ха го (2); 
ді дФ ді? добі ду? тд тд ба 
де  у- щільність середовища; с - питома теплоємність; 0 А- коефіцієнт 


теплопровідності; Т(р, ф.1) - температура середовища; і -час; т, - час релаксації. 

Математично задача визначення відносної температури циліндра (р, Ф, 1) 
складається в інтегруванні гіперболічних диференціальних рівнянь теплопровідності 
(2) в області Д, - Ко, фр є (рульря ) Фе (0,27), сє (0, го) , що, з урахуванням 
прийнятих допущень, запишеться у вигляді: 


90. | 026, 920. 41920. 100. 1 029, 
зад Я 


5 и ЕТ рт щу - 3 
о ор ай а ар рр рот| 
з початковими умовами 
28. (р,ф.б 
9. (р.р.0)- 0, спо) -о, 4) 
Ї 
граничними умовами 
в (ро) - Мр), 65 (і р.г)- У(ф), (6) 
умовами ідеального теплового контакту 
0 (р,.ф)- 65 (рі, фмг), (6) 
Я 20 (рьфаі) з, 90, (рі, фі) - (7), 
др др 
К Ко зі Р ж жі 5 - 
де р- в ;рог ГІ ро - КІ - радіус межі шарів; К, - внутрішній радіус циліндра; 
А коефіцієнт теплопровідності, У, - щільність, с, - питома теплоємність, а, - сви 
РУД 
г и ; з - 2 а; . Ми ; 
коефіцієнт | температуропровідності 5-го | шару циліндра; а; - - 2-;  8-1,2; 


(о) У (р) є С(Д). 
Тоді рішення крайової задачі  (3)-(7) 9.(р.р.1) є двічі неперервно 


диференційованим пор, Ф, і в області Р і неперервним на р |З, тобто 
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д.(р.фі)є С"|(р)П С(Р), а функції У/(ф) У(ф), 80, (р,ф.т) можуть бути розкладені в 
комплексний ряд Фур'є |41: 


8, (р, Ф.г) М ба, (р, 1) 


М) ро У М, у-ехрілю), (8) 
У(ф) У, 
де 
з (р,г) 1 ол 9, (р, ф.г) 
Ув. «бен УУ(ф) і-ехр(-іпруав, 
ло 
У, У(д) 


де 0, (ри) « 0) (руг)чні 00 р,ї), М, «Учні МО), му, ем ні МО), 
1- уявна одиниця. 

З огляду на те, що д.р.) функції дійсні, обмежимося надалі розглядом 
Є (р.г) для п-0, 1,2... тому що 8, ,, (р.г) То РР (р.г) будуть комплексно спряженими 


14). Підставляючи значення функцій з (8) у (3)-(7) у результаті одержимо систему 
диференціальних рівнянь: 


28 ддіту г 282. 0 зд | з/д 1980, п із 
ді п У5,п г 2 гУп ді шо: до? д 2 75,п 
р р ор 
з початковими умовами 
| 209(2,0 
9 (р.0)«0 ан а 110) 
граничною умовою 
00 (роугт)- МИ, | еф г)-У|), (1) 
умовами ідеального теплового контакту 
в (рг)є (ріг) (12) 
20 (р1 20 д1 
Й р) іі рт) з 


др Ор 
де 319 з-ат 92 з ат) т 22, то «Її; із1,2. 


Для розв'язання крайової задачі (9)-(13) побудуємо інтегральне перетворення: 


Нр НК 
й 1 (а 2) й 08 2) 
Ли) РУ р рмр з У, | 9 р(рідр (14) 
Ро 551 Р 5 
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пк 
"Б 2) ай якщо р є (ро Рі) 
1 


б "Ро 2) Р 


5 Мк 2 


2) 5 якщо ре (рі,р») 
02 


й 


Власні функції О, ем 
(24 


і власні значення знаходяться із ОЗв яЯзЗк задачі 
п 5 У 
5 


Штурма-Ліувілля: 


2 2 2 
і но "0,20 (15) 
ар" рр р' а 


и и 
"Р доро о о 16) 
711 2 


и и дО, р РІ | 00, Є РІ 
"РОЛИ й а, 2 2 (за) ЦТ) 
27 (27 р 


пк ИН 
ор о росту 
де дез 2) че ня | а 2) М а рову а 2) 
2 о 24 74 М 
Ник НК Ник НК На 
Чаз 2) - 3 , рр , 2) - И| , ре , 2) 
2 2 2 2 27 


129 


і (х) У, (х) - функції Бесселя 1" і 27? роду п"? порядку відповідно до (31. 


Власні значення й, Знаходяться із розв'язку трансцендентного рівняння: 
НК НК 
' , , 
ГУРА РІ Н РІ 
Гри 2 


1 
а р ЕЕ р 
24 2 
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де на 2) жи пк ЕЕ ре г я ре 2) са й 
2 2 2 2 2 2 М 


Формула оберненого перетворення має вигляд: 


"р)- У -Оиниат) ГАРНО; (18), 
пе 05 (шум) 


п 2 : 
де квадрат норми власної функції Ге МРУИ р) дорівнює: 


2 
- й, па |, 
Р) 75 282 2 0 
ОА Й кі 


- 19, й 
ца | 204 | тя | 
02 
а 
й, - й, і 
9, 7 А їз, 9 2.2 ад) 
- Р. 22 ав ШИ п - 22 
205 17 208 2 й 
В РР з Кр і Мак ПРярь м Р 
2 


Застосовуємо до системи диференціальних рівнянь (9) інтегральне перетворення 
(14) і, враховуючи позначення (1), в результаті одержуємо систему звичайних 
диференціальних рівнянь: 


а (т,) 24(ї) | 
ч ЗСУ І | І о нь Об) -м2810 з (19) 


аб 
І аї ад | 


де ОЇ - то бай о з 2) УФ 
1 


74 


з початковими умовами 


а (б) 
(ши) є 0, сб, люде о ооваб 020) 


Застосовуємо до системи диференціальних рівнянь (19) з початковими умовами (20) 
інтегральне перетворення Лапласа |41: 


ор 
75) - Гс)е"?? г 
0 
У результаті одержуємо систему алгебраїчних рівнянь відносно 81) 
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збо пробій | з от 


де 1-1,2; чик З зай - 


Розв'язавши систему рівнянь (21), одержуємо: 


сі) 2 ін б/ть ;) 
ба но РУС у (- 1" ат Ок із.) 


а-с1,2) (22), 
(за здо на? п Мічзт, є 


де пк - Мп 


Застосовуючи до зображення функцій (22) формули оберненого перетворення 
Лапласа, одержуємо оригінали функцій: 


90) б. г)з Хол, ус) (з; )|оеуз; «тупі 025; а ам -фе, 5і ЇЇ 
ен Хоо УВУ Мало зано Мото 
ри ) Р (23) 
в) сао МО, Мен тні В Ми ун 
у чі іа. (з бе), У ре,з; таті) - 0) (з; ) рат (дез нї) 


ром) 04) 


-1 
0.55; ск 


, а значення 5; для 1-1,4,3,4 визначаються за 


де і )- (туз; зар -(гуап) 


формулами 


(гуаті - туя а з туапі)? 474 (гуаті . туя а - туапі)" 474 
ЧР ог б ог | 
ТД т 


Таким чином, з урахуванням формул обернених перетворень (8) 1 (18), одержуємо 
температурне поле кусково-однорідного кругового циліндра в напрямку полярного 
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радіуса, який обертається з постійною кутовою швидкістю () навколо осі О7, з 
урахуванням кінцевої швидкості поширення тепла: 


НК 
бан вен Й "Р 
Жоа)ч 5158 или) | 2 
кози о 2) 


С 


«ехр(іпф), 


де значення 8. к .1) 180)шид .) визначаються за формулами (23),(24). 


Висновки 

У статті, за допомогою розробленого нового інтегрального перетворення, 
знайдено температурне поле порожнього кусково-однорідного кругового циліндра в 
напрямку полярного радіуса, який обертається з постійною кутовою швидкістю 
навколо осі 07, з урахуванням кінцевої швидкості поширення тепла, у вигляді збіжних 
ортогональних рядів по функціях Бесселя і Фур'є. Знайдено аналітичний розв'язок 
узагальненої крайової задачі теплообміну порожнього циліндра, який обертається, з 
урахуванням скінченності величини швидкості поширення тепла, і може знайти 
застосування при модулюванні температурних полів, які виникають у багатьох 
технічних системах (у супутниках, прокатних валках, турбінах 1 т.1.). 
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КЕ5ОМЕ 

М.С. Вегапук 

Машетабісаї піодеї апа тепод ої 5о0Їийоп уїі Ше пеуу іпіесга! кгап5їогтабйоп 
ої а гепегаїгед Боппдагу ргобіет ріесеууї5е ппіїогт Неаї ігап5їег суппаєг 

Тріз5 рарег деаї5 мйіїб Ше саїсціайоп ої ип5іеаду іеппрегаїиге Пе егаріу суппаєвг 
роПом сігсиіаг суппаєг оціег гайіця К їп а суппагіса! соогадїпа/е зузіет, ріесем/т5е ппіїогт 
іп Фе дігесйоп ої Фе роіаг гадїця г, місії гокагез ас а соп5іапі апеціаг уеЇосііу ахі5 О2, 
сакіпє їпіо ассойпі (Пе Ппіїе 5рееад ої ргорабайоп ої реак. Трегтаї ргорегіїе5 аге іп еаср 
Іауег іпдерепдепі ої (еппрегаките сопаїйоп5 1деаї (ШФегтаї сопіасі Бесууееп Ше Іаусг5 ап 
іпіегпа! Беас 5о0игсе ауаїйабіе. Ас Фе іпійа! йте Ше суппаег іетрегаїиге сопзіапі, апа Фе 
оціег апа іппег 5игіасез ої Фе суппдаєвг іетрегагиге 15 Кпомуп апа дое5 пої Ферепа оп йте. 

Тре тафетайса! пиодеї ої (еппрегагиге дії5пібийоп іп а ріесеу/ї5е ПБотоєепесиця 
суппдег, урісі го(аге8 ай а соп5(апі апеціаг уеЇосіїу ахі5 О7 (акіпє їпіо ассоцпі Фе Ппіїе 
зреед ої ргорагацоп ої Пеає іп Ше їоги ої птааїбетайса! рпузіс5 Бойпдагу ргобіет Їога 
зубіет ої Пурегропс даїегепиа! гіупуапуап ШФегта! сопдйсйуїсу ул бойпдаагу сопдаїйоп5 
ої Фе Нг5і Кіпа. То 50Їме Ше Бойпдаагу ргобіеті гесеїуед Ше геадиігед (етрегаите Пе!д 15 
рге5епіса а5 а соппріех Коигіег 5егіе5. Реуеіорей пему іп(еєта! (гап5огтайоп Їог ріесем/і5е 
ротоєепеойя 5расе м/р у/бісп Їоцпа (Ше (еттрегаїиге Пеід ої сопйпиоц5 ріесему/іве 
ротоєепеоця сігсціаг суппаєг іп а сопуегеепі 5егіе5 іп 5ресіа! огірборбопа Гипспоп5 апа 
Еоигіег Везееі. 

ЕБоипа апаїуйасаї 50Їайоп ої Фе гепегайгед Бойпдаагу ргобіета ої реа ігап5їег суппаєвг 
аг гокае5, бімеп Ппіге уеіосіку ої ргорабацоп ої Беає сап Бе ц5ед іп Ше подніацоп ої 
іетрегаїиге Пе!45 агізіпе їп плапу КесПпіса! 5у5(етзя (п 5ахеШтез, го8, кагбіпез, еїс.). 
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